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Weyl semimetals typically appear in systems in which either time-reversal (T) or inversion ( P ) 
symmetry is broken. Here we show that in the presence of gauge potentials these topological states 
of matter can also arise in fermionic lattices preserving both T and P. We analyze in detail the 
case of a cubic lattice model with 7r-fluxes, discussing the role of gauge symmetries in the formation 
of Weyl points and the difference between the physical and the canonical T and P symmetries. 
We examine the robustness of this PT invariant Weyl semimetal phase against perturbations that 
remove the chiral sublattice symmetries and we discuss further generalizations. Finally, motivated 
by advances in ultracold atom experiments and by the possibility of using synthetic magnetic fields, 
we study the effect of random perturbations of the magnetic fluxes, which can be compared to a 
local disorder in realistic scenarios. 


I. INTRODUCTION 

In recent years it has been understood that, besides 
gapped phases of matter 1 ®, gapless systems may also 
present topological properties. The main examples are 
topological metals and semimetals, whose band structure 
reproduces the main features of topological insulators 
and superconductors at their critical points. This is the 
case of the so-called Weyl semimetald3®, whose recent 
experimental realizations include tantalum arsenidd-^®, 
niobium arsenide®! and photonic crystal^®. 

A crucial ingredient in obtaining Weyl semimetals 
is the br eaking of inversion (P) or time-reversal (T) 
symmetry ! 5 I 1 A \15l j n the absence of at least one of 
these symmetries it is possible to obtain zero-energy 
bulk modes with linear dispersion (Weyl nodes), which 
constitute monopoles for the Berry connection and are 
therefore protected against generic perturbation^®. The 
nodes lead to the appearance of zero-energy surface states 
(Fermi arcs), connecting the projections of Weyl points 
with opposite monopole charges on the surface Brillouin 
zon^ (see Fig. [l]). 

If the Hamiltonian of the system is invariant under 
both T and P, starting from a Weyl node at momentum 
ko and applying T one obtains a node with the same 
chirality at — fco- The subsequent application of P creates 
a new Weyl point with opposite chirality in the original 
point k 0 . As such, one cannot isolate single monopoles 
of the Berry phase® (see ref. [TTJfor a similar analysis on 
Dirac semimetals). This mechanism, however, relies on 
the fact that both T and P commute with the translation 
operator. 

Here we show that a different scenario can occur when 
gauge symmetries are present. In this case, the previous 
argument no longer holds because the Weyl points linked 
by T and P do not appear in —ko and fco- In the pres¬ 
ence of a gauge potential A , we must distinguish between 


a momentum k and a physical momentum k+A: only the 
latter is inverted by T and P. Therefore, one must distin¬ 
guish “physical” T and P symmetries, which are realized 
with the addition of space-dependent gauge transforma¬ 
tions, and ’’canonical” T and V symmetries, expressed in 
terms of translationally invariant operators, mapping the 
Hamiltonian to a different, but gauge-equivalent, form. 
This distinction is similar to the one between magnetic 
and canonical translation operators. 

This observation is motivated by the recent experimen¬ 
tal advances in ultracold atom setups. Several experi¬ 
ments showed that it is possible to obtain large artificial 
magnetic fluxes in optical lattices through laser-assisted 
tunnelingf®!®. This offers new possibilities to reach a 
Weyl semimetal regime based on cubic lattices subject 

tO 7T- fluxepHU, 

not breaking the physical T and P sym¬ 
metries. The experimental observation of Weyl points 
in these setups can be achieved through band mapping 
techniques, based, for example, on Bragg spectroscopy®. 
An alternative detection procedure relies on the Berry 
curvature of the energy bands, recently obtained in two- 
dimensional setups through interferometric techniques 26 
and measurements of the band population^®. These 
methods can generalize the results obtained for two- 
dimensional Dirac pointd® to three space dimensions. 

Our main conclusion is that, in the presence of gauge 
symmetries, one needs to break VT to have Weyl points, 
but PT can be preserved. To derive our results, we ex¬ 
amine first the role of gauge symmetry in the cubic lat¬ 
tice model with magnetic 7r-fluxes. It describes a three- 
dimensional generalization of the experimental results 
obtained in Refs. ITSl - Ef)! and can be engineered by op¬ 
tical toolj23 29 3D . This model illustrates the simplest 
occurrence of Weyl semimetals in a PT invariant sys¬ 
tem with gauge invariance. The robustness of this phase 
against the breaking of the chiral sublattice symmetry is 
shown. Later on, we study the effect of perturbations 
in the fluxes, which may characterize the physical real- 
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ization of this Weyl phase. Such perturbations, despite 
breaking of PT, do not spoil the Weyl physics for realis¬ 
tic system sizes and sufficiently small random variations. 
Finally we investigate the effects of a parabolic trapping 
potential on the Weyl cones and on the related Fermi 
arcs, showing that they survive, even if deformed, in this 
condition. 


II. CUBIC LATTICE MODEL WITH n FLUX 

Consider a tight-binding model of spinless fermions 
on a 7r-flux cubic lattice (we set the lattice constant 
to a = 1). Given a specific gauge, such a system 
breaks either the canonical or 7 HH symmetry, de¬ 
spite its physics being P and T invariant. We work in 
the Hasegawa gauge 2 -^, with a vector potential: 

A = ir(z-y,y-z,0), (1) 

corresponding to a magnetic field B = n (1,1,1). The 
tight-binding Hamiltonian of the system under the gauge 
potential Q reads 

H({6}) = -w Y, c\ +f ie ^c P + H.c., (2) 

r,j 

where c 1 and c are creation and annihilation operators 
for fermions on the lattice, u> is the hopping strength, 
j = x,y,z, and Qj{r) = J^ 3+a Aj (r) drj. Eq. (jTJ) leads 
to: 


0 X = 7T {z - y) , 6 y = tv (y - z) + 7t/2 , 9 Z = 0 . (3) 

This gauge potential is V invariant: an inversion trans¬ 
formation centered in a lattice site, 

x->-x, y->-y, z-t-z, (4) 

leaves all phases 9j unchanged, since e l ' K ^ z ~ v ' > = 
e -m(z-y) Therefore the physical and canonical inver¬ 
sion symmetries are equivalent with this gauge choice, 
P = V, and H = Pi HP. 

The situation is different for T. We define it as com¬ 
plex conjugation of the hopping amplitudes, such that 
TH({9})Tl = H({—9})-, this relation implies that T is 
broken with our gauge choice, since e tey< ' v,z ' > 7 ^ e _z T<T 2 )_ 
The system is however T invariant. Indeed, through a 
gauge transformation U y c? = e lny c?, H({—9}) is equiva¬ 
lent to H({9 }): 

TH{{6})Tl =UlH{{6})U v . (5) 

We conclude that the physical time-reversal transfor¬ 
mation, which leaves the Hamiltonian in Eq. ([ 2 ]) un¬ 
changed, takes the form T = U y T , requiring a space- 
dependent gauge transformation to be implemented. In 
contrast, the space inversion V = P has a canonical, 
translational-invariant form with our gauge choice. A 



FIG. 1. Left: Effect of time-reversal and inversion symme¬ 
tries on the Weyl cones of the Hamiltonian ©• A Weyl cone 
is mapped onto itself when applying the canonical VT sym¬ 
metry. However, due to the gauge transformation U y , the 
physical PT symmetry relates Weyl cones at different mo¬ 
menta. Right: Bandstructure in a diagonal slab geometry, 
with hard-wall boundaries in the x — y direction and a thick¬ 
ness of 80 sites. The color scale indicates the amplitude of the 
eigenstates in the first and last eight sites from the boundaries. 
Weyl cones of opposite chirality are connected by zero-energy 
Fermi arcs (in red). 


different gauge choice can lead to the opposite situation 
where T = T and V 7 ^ P (see ref. [231 and Appendix [A]) . 

Our main point is that, no matter the gauge choice, 
the system is not PT invariant, which is necessary to 
have Weyl points, but actually it is invariant under the 
physical PT symmetry, even though a space-dependent 
gauge transformation is required to make this invariance 
explicit. 

The effect of U y is to flip the signs in the tunneling 
amplitudes of Eq. 0 along the y direction, and, given 
its local nature, can be applied independently on the 
boundary conditions of the system. This is a general 
feature of every gauge transformation. For our particu¬ 
lar gauge choice in Eq. 0 the physical P and T sym¬ 
metries commute; the same happens with every gauge 
choice because, given a unitary transformation W, one 
has [WtPW, WtTW] = W t [P, T\W = 0 . 

Following refs. [21! and [51;, we rewrite Eq. ([ 2 ]) in mo¬ 
mentum space: 

H = —2 u) (t x cos k x — t v cos k v + t z cos k z ), ( 6 ) 

where t* are Pauli matrices and k belongs to the mag¬ 
netic Brillouin zone (BZ). In the gauge ([I]), we take 
k x , k y £ [—7r,7r] and k z £ [0,7r], and identify ( k x ,k y ,k z ) 
with ( k x , k v + 7T, k z + 7rp2l. The eigenstates of t x corre¬ 
spond to the even and odd sublattices in the y — z plane. 

The energy bands of the system touch at the Weyl 
nodes k^'^ 1 = (±7r/2, ±7 t/2, 7t/2 ), having a chirality 
x(^o = ’ ± ) = Uj= x , y , z ^ cos%|£ = g±,±, such that the Weyl 

points at k^ ,ziz are monopoles of charge x f° r the Berry 
connection. 

In Eq. ([ 6 ]), T does not take a canonical form, al¬ 
lowing the formation of Weyl points. A Weyl point 
at is mapped by T = U y T into a point k' 0 = 
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(—7r/2, — 7r/2 + 7r, — 7 t/ 2) with the same chirality, corre¬ 
sponding to kf'~. V transforms fc ( 7~ into a node with 
opposite chirality at — kff~ = k^~. Therefore, the over¬ 
all transformation VTU y relates inequivalent Weyl points 
with opposite chirality at different momenta, and 

as shown in Fig. [l] 

The PT invariant Weyl semimetal is a universal phase 
emerging from different lattice geometries with magnetic 
fluxes. In Appendix [B] we present a family of Hamiltoni¬ 
ans interpolating between the 7r-flux cubic lattice model 
and one defined by a stack of horizontal brickwall lat¬ 
tices, connected by vertical square plaquettes with n- 
fluxes. The models in this family are PT invariant and 
display Weyl points. The definition of their physical sym¬ 
metries relies, in general, on gauge transformations, and 
their effect on the Weyl cones can be analyzed following 
the general ideas presented above, leading to conclusions 
analogous to the ones depicted in the scheme in the left 
panel of Fig. & (see Appendix IB]). 

Weyl nodes have a topological origin, since they are 
monopoles of the Berry curvature for the energy bands. 
This feature induces stability against local noised How¬ 
ever, perturbations which couple two Weyl points with 
different momentum and opposite chirality cause them 
to annihilate pairwise. In the following we discuss the 
stability of the PT invariant Weyl phase against pertur¬ 
bations breaking the chiral sublattice symmetry or con¬ 
sisting of small fluctuations of the fluxes in finite-size 
systems. 


Breaking of the sublattice symmetry 

The Hamiltonian in Eq. (pi) is characterized by a chiral 
sublattice symmetry c ? —> ("7-1 ) x+v+z cp mapping H into 
—FT, possibly broken by the introduction of staggering 
perturbations or onsite disorder. Below we show that 
this symmetry does not play any fundamental role in 
protecting the Weyl points: PT invariant Weyl phases 
can exist also in systems without it. 

We consider two staggering on-site potentials: Vtj = 
V (—l) ri_r J and V xyz = V {—l) x+v+z , both breaking the 
chiral symmetry. In their presence, the BZ is further 
halved along two directions. Correspondingly, the spec¬ 
trum of Eq. (p]) acquires additional degeneracies and 
Weyl points with different chiralities overlap in momen¬ 
tum space. A moderate perturbation 14, does not open a 
gap: it removes one of the additional degeneracies and it 
gaps only half of the Weyl points, leaving a pair of (degen¬ 
erate) Weyl cones. Therefore the system remains in the 
PT invariant Weyl semimetal phase, despite the break¬ 
ing of the sublattice symmetry. The potential V xyz , in¬ 
stead, immediately opens a bulk gap^ by coupling Weyl 
nodes with different chiralities. However, this coupling 
is avoided when the Weyl points are displaced in mo¬ 
mentum space by a term V % j, stabilizing the Weyl phase 
against V xyz as long as |I4, | > |I4y2|- We conclude that, 


even in the absence of a chiral sublattice symmetry, there 
exists in general an extended Weyl semimetal phase in¬ 
variant under the physical P and T symmetries, and sta¬ 
ble against onsite perturbations not breaking them. 


III. PERTURBATION OF THE FLUXES 


A deviation of the fluxes from 7 r can couple separated 
points in the BZ; therefore it is important to study its 
effect on the Weyl phase, crucially evaluating the range 
of errors in the fluxes still enabling the observation of the 
Weyl cones. We consider then the vector potential: 


A = (/3z — 72/,7T2/ - az, 0), 


(7) 


corresponding to B = (a, 4 , 7 ). This configuration im¬ 
plies in general a breaking of T and P, driving the system 
from the PT invariant point to one with no symmetries. 
Indeed, both T and P invert all the magnetic fluxes, thus 
mapping the model into a gauge inequivalent system, un¬ 
less a = /? = 7 = 7T. 

We consider flux configurations 4 = tt(1 — l/g^), with 
(f> = a,/ 3,7 and q$ being a triplet of integers, leading to 
an increase in the magnetic unit cell, now whith a vol¬ 
ume given by the leas t com mon multiple (LCM) of qw. 
V uc = LCM(q a ,q/3,q~ / '^2fM. We find that in the ther¬ 
modynamic limit flux perturbations can have a drastic 
effect on the low-energy physics of t he sy stem, also due 
to the fractal nature of its spectrunPSES For instance, 
when q a = qp = q 1 > 2 the system is in a metallic 
state at half filling, with a two-dimensional Fermi sur¬ 
face. A different case occurs for anisotropic configura¬ 
tions B = (7r(l — l/q),7r,7r): for even q the spectrum 
shows Weyl cones, whose density in the magnetic Bril- 
louin zone increases with q. This feature requires larger 
and larger system sizes to be seen for decreasing per¬ 
turbations l/g. For example, in Fig. [ 2 ] we show the 
spectrum for a = 37r/4 and 4 = 7 = 7 r. Interestingly, 
the physics of these Weyl points appears only at rel¬ 
atively low energies, characterized by a threshold ew, 
then a sufficiently high energy resolution is required to 
detect the Weyl cones. For energies between ew and the 
bandwidth, instead, the system is well approximated by 
a metallic behavior. For these anisotropic flux configu¬ 
rations the chiral symmetry constraints the Weyl points 
to be at zero energy, even in the presence of a simultane¬ 
ous breaking of T and P, differently from the standard 
behavior in other setups 37 -*^. Finally, for odd q the low 
energy physics is described by a node-line semimetaP^ 
with a one-dimensional Fermi surface, instead of a Weyl 
semimetal. 

In ultracold atom experiments, lattice sizes are typi¬ 
cally limited to tens of sites. Therefore, for small pertur¬ 
bations from a = /3 = 7 = 7 r, the available system sizes 
are much smaller than the magnetic unit cell, V <C V uc . 
Below we show that, in this limit, the Weyl physics is ro¬ 
bust to small variations of the fluxes. To examine these 
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FIG. 2. Spectrum for B = (3/4,1, l)7r with k z = 0, displaying 
Weyl cones at zero energy. The inset shows the detail of a 
Weyl cone: the topological semimetal physics emerges clearly 
below an energy scale ew , whereas above it the system has a 
metallic behavior. 


perturbations, we compute the density of states (DOS) 
of a cubic lattice with a finite linear dimension L, using a 
kernel polynomial approximation (see Appendix [C]). The 
DOS is averaged over random flux configurations where 
a, /3, 7 are independent variables of the form (j> = 7 r + <5</>, 
and S(f> is drawn randomly from the uniform distribution 
[— {/, U], with U being the maximum perturbation. In 
this way, both the strength and the direction of B are 
varied. The presence of a Weyl semimetal phase can be 
determined by examining the energy dependence of the 
average DOS, taking the form p(E ) ~ E 2 for linearly 
dispersing states. 

We investigated different values of U up to 0.05, for 
different system sizes up to L = 160. Our results for the 
DOS share similarities with the physics of Weyl semimet¬ 
als in the presence of onsite disorder discussed h j 40 * 41 ( 
For relatively small values of U, the zero-energy den¬ 
sity p(E = 0) is different from zero and, in agreement 
with recent numerical results®, two energy thresholds 
emerge: such that the average DOS is approximately 

constant for \E\ < ei, and ew such that p(E) oc E 2 for 
£i < \E\ < ew, as expected for Weyl points. 

The threshold e\ sig nals in general the appearanc e of a 
diffusive behavio j 42 * 43 ! at low energies. Recent works®® 
suggest that this constant DOS can be associated with 
rare power-law localized states, although for flux per¬ 
turbations its origin may be different. Independently 
on their nature, these states may hinder the observation 
of Weyl quasiparticles: the visibility of the Weyl phase 
depends in general on the width of the energy window 
ew ~ £i, which must be significantly large. 

The average DOS for L = 160 is plotted in Fig. [3] For 
U = 0.01, we find e\ ~ 0.3w and ew ~ w. Close to the 
Fermi level the DOS takes an energy independent value 
p(E/u> = 0) « 10~ 3 , whereas at intermediate energies 
we recover the scaling associated with a Weyl semimetal 
phase. Using the fit p(E) — p( 0) = a\E\ b for U = 0.01 in¬ 
side the interval [e\,ew] (we choose 0.52 < \E\/uj < 0.75) 
we obtain a ~ 0.029 and b ~ 2.12, with a reduced 
X 2 — 1.02. For smaller system sizes the qualitative 



FIG. 3. Double logarithmic plot of the density of states, as 
a function of energy for different values of U. The DOS is 
computed for a cube with linear size of 160 sites, and is aver¬ 
aged over ~ 150 — 200 random flux realizations. For clarity, 
the green and blue curves are shifted vertically by multiplying 
with a constant factor (shown on the plot). The DOS close 
to the Fermi level is constant, indicating a diffusive metal 
behavior. In the intermediate energy range and for small val¬ 
ues of U we observe an energy scaling characteristic of Weyl 
semimetals (dashed lines), p(E) — p(0) = a\E\ b with b « 2. 
Larger values of U lead to an increase in the energy range as¬ 
sociated with a diffusive metal, and a reduction of the region 
associated with a topological semimetal. 


behavior is very similar, even though the average DOS 
shows pronounced finite size effects, leading to a larger 
value of x' 2 . As U is progressively increased, the energy 
ei associated with diffusive metal behavior also increases, 
while the topological semimetal region between e\ and ew 
becomes smaller, eventually disappearing (as shown for 
U = 0.05 in Fig. [3]). 

Our results indicate that the Weyl semimetal behavior 
is observable in ultracold atom implementations of the 
7r-flux cubic lattice model, given variations in the flux of 
1% — 2% for L = 160 (a realistic estimate for setups of 
this kincpHlol) 

and generally 5cj) < 7 t/L. 


IV. EFFECT OF A CONFINING POTENTIAL 

Among the perturbations relevant for ultracold atom 
setups, which may have an impact on observing Weyl 
points and Fermi arcs, the presence of trapping potentials 
deserves particular attention. 

Experimental setups in optical lattices usually rely on 
harmonic traps or hard-wall potentials obtain ed th rough 
optical box trapJ® or light-intensity maskiP^HZ], and, 
typically, they reach a system size of the order L ss 100. 
We examine here the fate of the zero-energy mode^ in 
the presence of both these potentials. 

Since the presence of a boundary along the y axis does 
not break the gauge invariance of Eq. (2) of the main 
text, it does not spoil the distinction between canonical 
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FIG. 4. States with energy |E| < 0.02w in the presence of 
a harmonic trap in the surface BZ. The color depicts their 
localization close to the hard-wall boundary. Here L = 50, 
y = ui/2 and is chosen to have the potential energy rang¬ 
ing from —uj/2 to +w/2 from the center towards the surface. 
The Weyl points become extended regions of zero-energy bulk 
states, connected by curved Fermi arcs. 

and physical symmetries, relying on the U y gauge trans¬ 
formation. One consequence is that the same boundary 
can still host Fermi arcs, as in different geometries not in¬ 
volving the y axis. Keeping in mind this observation, we 
consider a diagonal slab geometry with hard wall bound¬ 
aries located at \x — y\ = L , and a harmonic trapping 
potential depending only on the distance from the sur¬ 
face. The single-particle Hamiltonian reads 

H con{ = H + y fl 2 (x - y) 2 - y, (8) 

with H as in Eq. (2). For a small chemical potential, 
0 < y < 2 w, a metallic phase appears at x = y, with 
the second energy band partially filled. In a local ap¬ 
proximation, the zero-energy Fermi surface corresponds 
to contours around the Weyl points. Moving towards 
the surface, the harmonic trap reduces the local chemi¬ 
cal potential. Therefore, the Fermi surface shrinks into 
the Weyl points, which are reached at a distance 1Z such 
that m VL 2 TZ 2 = 2y. Further increasing \x — y\, the sys¬ 
tem is driven again in a metallic phase, and, reaching the 
surface at \x — y\ = L , zero-energy surface modes appear. 

Figure [4] shows the states close to zero energy on the 
surface BZ in a diagonal slab geometry. As expected, the 
zero-energy states of the full system include two regions 
of bulk states surrounding the projections of the Weyl 
points. Surface modes interpolate, with a curved shape, 
between the Fermi surfaces around the Weyl points. 
These Fermi arcs can be experimentally probed through 
a band mapping based on Bragg spectroscopjES. 


V. OUTLOOK 

In fermionic systems with gauge invariance the phys¬ 
ical time-reversal and inversion symmetries must be de¬ 
fined up to local gauge transformations. Such transfor¬ 
mations have the non-trivial effect of translating states in 
momentum space, allowing to avoid the usual doubling 
of the zero-energy Weyl points present in systems with 


canonical time-reversal or inversion symmetries. There¬ 
fore, local gauge symmetries can be exploited to define 
a new kind of PT invariant Weyl semimetal. We dis¬ 
cussed this phase for the simplest case of a cubic lattice 
model of spinless fermions with 7r-magnetic fluxes, show¬ 
ing that its Hamiltonian is invariant under a transforma¬ 
tion VTU = PT, despite the presence of Weyl points in 
its spectrum. We emphasize, however, that the existence 
of the PT invariant Weyl phase is not limited to spinless 
fermions with a physical time-reversal symmetry obey¬ 
ing T 2 = 1, but it can be easily extended also to spinful 
models with T 2 = — 1 (see Appendix |d| . 

Furthemore, motivated by the advances in ultracold 
atom experiments where synthetic gauge fields may be 
added and tuned, we discussed the robustness of this 
topological semimetal phase against the breaking of the 
chiral lattice symmetry, the fluctuations of the flux for 
realistic system sizes, and the introduction of a parabolic 
trapping potential. 

Our results suggest that the role of gauge symme¬ 
tries in topological phases of matter have been so far 
underestimated. The tenfold classification of topologi¬ 
cal insulators and superconductors 1 " 3 is based on non- 
spatial (time-reversal, particle-hole and chiral) symme¬ 
tries, which are resilient against most sources of disor¬ 
der. This classification has been extended to spatial 
symmetries®, including poi nt gro up symmetric ^® 0 and 
their magnetic counterparts- 152 . However, gauge sym¬ 
metries cannot be reduced to them and a systematic 
study of their effects is still missing. Also for gapless 
systems a complete classification of topological phases 
exists based on non-spatial symmetrie d 3 ^ 53 ^ and only a 
few generalizations have been considered so fai^. 
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Appendix A: Alternative gauge choice 

The gauge choice used in the main text is invariant 
not only under the space-inversion symmetry centered on 
the lattice sites, but also on space inversion symmetries 
centered on the lattice cubes, on the plaquettes in the yz 
planes and on the links in the x direction. Different gauge 
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choices break some of these symmetries. For example in 
Ref. [24l Dubcek et al. use: 

Ox = -7T (x + y) , 9 V = 0, 6 z =n{x + y). (Al) 

In this case H ({0}) = H ({—0}), such that T = T. In 
contrast, the effect of the space inversion centered on the 
links in the x direction is: 

e"--»•-e w * , e w ‘(A2) 

such that the Hamiltonian in momentum space is not 
invariant under the inversion k x —> — 

T(fc) = — 2 w [(Tj; cos(k y ) + <j y sin(fc x ) + cr z cos(fc z )]. 

. (A3) 

In this case a gauge transformation U xz = e l7r ( x+z ) 
is required to define the physical inversion symmetry: 
P = VU XZ . Thus, the PT transformation takes the form 
PTU XZ , which confirms that the system is PTinvariant 
but not PTinvariant. 


Appendix B: A generalized model for the PT 
invariant Weyl semimetal phase 


The cubic lattice model that we presented in the main 
text is only a specific example of a PT invariant Weyl 
semimetal phase. In this section we discuss first an al¬ 
ternative model, based on a stack of brick-wall lattices, 
that hosts the same gapless topological (Weyl) phase pre¬ 
serving the physical PT invariance; then we show that it 
is possible to interpolate continuously between the cubic 
lattice model analyzed in the main text and this alterna¬ 
tive lattice. The alternative model is based on a three- 
dimensional lattice defined by brick-wall layers piled such 
that all the lattice sites in different layers, labeled by z, 
overlap when projected on the xy plane. Therefore all the 
plaquettes oriented in the xy planes are rectangular with 
dimension 2 a x a, whereas all the plaquettes in the xz 
or yz orientations are square plaquettes. We introduce a 
magnetic 7 r-flux in only the vertical plaquettes. In this 
way, such a three-dimensional lattice can be obtained by 
suppressing half of the tunnelings along the y direction 
in the cubic lattice model with 7r-fluxes in all the direc¬ 
tions. The horizontal planes xy, indeed, are constituted 
only by rectangular plaquettes in which the magnetic flux 
is doubled to 27 t and, therefore, has no physical effect. 

The real space Hamiltonian results: 


Hbw = -W [e w » (i H +c °r + H.C.' 


r : x—y=eve n 


'E + e ^ (f H+* W + H.c 


(Bl) 


where we can choose a gauge such that: 

Q z (r) = n(x - y ), 0 X = 0 y = 0 . 


(B2) 


Considering a two-site unit cell given by the sublattices in 
the horizontal planes, and labeling it by the pseudospin 
t z = ±1, we obtain the following Hamiltonian in momen¬ 
tum space: 

Hbw 

-= — 2 r z cos(fc z )— 2t x cos (k x )—r x cos(k y )+T y sin (k y ). 

UJ 

(B3) 

The related spectrum results: 


E = ±w4cos 2 k z + (2 cos k x + cos k y ) 2 + sin 2 k y (B4) 

and it is characterized by four Weyl points in = 
(±27t/3, 0, ± 7 t/ 2 ) with a choice of the BZ such that 
k x ,k z £ [— 7r,7r) and k y £ [— 7t/2, 7 t/ 2 ). 

The previous Hamiltonian is trivially invariant under 
the canonical time-reversal symmetry T, since the Hamil¬ 
tonian (Bl I has only real hopping amplitudes. Therefore 


T = T. Besides, the lattice is invariant under a space 
inversion centered in one of the rectangular plaquettes, 
such that V is described by: 


y -a i - y, 


z —z. 


(B5) 


However, our gauge choice is such that 9 Z is not invariant 
under this canonical transformation and e l6 - —e l6 -. 
Therefore the physical space inversion must be defined by 
P = PU Z , where the gauge transformation U z c? = e I7rz cp 
enters. 

It is easy to see that T maps fco’ b whereas 

P = PU Z gives kQ ,b —> kQ a,b . Therefore also in this case, 
despite the PT invariance, Weyl points with different chi¬ 
ralities have different positions in momentum space, thus 
ensuring the possibility of the Weyl phase. 

It is possible to interpolate between the brick wall lat¬ 


tice model in Eq. (Bl) and the cubic lattice model with 


7 r-fluxes discussed in the main text. This is achieved by 
adding the missing hopping terms along the odd links in 
the y directions of the first Hamiltonian. More formally, 
we define: 


H(a) = Hbw — Q-co 5 l 

r s.t. x—y —odd 


e ie ^4 + .cr + H.c. 


(B 6 ) 


We complete the definition of the phases in Eq. (B2) 
with: 


d y (r) = n(x - y). (B7) 

With this choice of the gauge, H(a) describes again a cu¬ 
bic lattice model with 7r-fluxes in all the plaquettes, as in 
the main text. It interpolates between H(a = 0) = Hbw 
and the Hamiltonian H(a = 1) which is analogous, up 
a rotation, to the cubic lattice model with the hoppings 
in Eq. ©• For each value of a the canonical T is 
preserved, whereas the physical inversion symmetry P 
coincides with the one specified above for Hbw- 

The position along k x of the Weyl points varies with 
a, = (±q(a), 0, ±7t/ 2), but the relation P|fcQ ,h ) = 
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|fc-“’ b > remains valid for each value of a € [ 0 , 1 ], such 
that all the Hamiltonians obtained in Eq. (B 6 ) describe 
a PT invariant Weyl semimetal. 


Appendix C: Kernel polynomial method 

In Fig. 3 of the main text, the density of states p{E) is 
computed using the kernel polynomial method®. In the 
following we summarize the basic steps of this approxi¬ 
mation. For a Hamiltonian H with eigenvalues E n , the 
density of states is given by 


the density of states plotted in Fig. 3 we have used the 
so-called Jackson kernel, corresponding to 


(N — n + 1) cos 


9n — 


7m 

N + 1 


nn 


sin 


cot 


iV + 1 N + l 


N + 1 


(C 8 ) 


and computed the first N = 2048 terms of Eq. (C7|. 

Finally, to determine the coefficients (C 6 ) we have used 
a stochastic evaluation of the trace 


Pn = Tr T n (H) 




(C9) 


p(E) = Y,S{E-E n ). 


(Cl) 


We will expand Eq. (Cl) in a series of Chebyshev poly¬ 


nomials, T n (x) = cos[narccos(x)], which obey the recur¬ 
sion relations 


T n+1 {x) = 2xT n (x) - T n _i{x), (C2) 


with Tq{x) = 1 and T\(x) = x. Since the polynomi¬ 
als are defined only in the interval [— 1 , 1 ], the first step 
is to rescale the Hamiltonian such that its spectrum is 
contained this interval: 


H = 


H-b 

a 


such that the energy becomes 


E = 


E-b 

a 


(C3) 


(C4) 


Setting uj = 1 in the Hamiltonian (2) of the main text, we 
choose the values 6 = 0 and a = 8 , for which the rescaled 
energies E £ [—1,1]- Following Ref. [56, the rescaled 
density of states, p{E) = ~ E n ), can be written 

as an infinite series 


p(E) 


1 

7T \/1 — E 2 


oo 

P o + 2 p n T n (E) , 

71 — 1 


with the expansion coefficients given by 
p n = Tr T n (H), 


(C5) 


(C 6 ) 


where Tr denotes the trace. 

For our numerical purposes, we truncate the expansion 


of Eq. (C5| by keeping only the first N coefficients: 


p{E)^ 


[Vi - E 2 


Po 


N—l 

2 El p n T n {E ) 


(C7) 


This approximation leads to fluctuations in the density of 
states, also known as Gibbs oscillations, which we damp 
by modifying the expansion coefficients p n —> g n p n ■ For 


where |r) are random states with entries drawn from the 
Gaussian distribution with zero mean and unit variance 
A/"(0,1). Each coefficient was computed using R = 10 
random vectors. 


Appendix D: Spin 1/2 PT-invariant Weyl semimetal 

It is possible to obtain a PT-invariant Weyl semimetal 
phase also in spinful models characterized by T 2 = — 1, 
differently from the spinless case of Eq. ([2]) for which 
T 2 = +1. To investigate this scenario, let us study a toy 
model of spin 1 /2 fermions on the cubic lattice in which 
both the spin species are subject to 7r-fluxes in all the 
plaquettes and they are mixed by non-trivial tunneling 
operators. We consider the Hamiltonian: 

H 2 ({9}) = ~u E cl +l O ss> eie ^ Cr ^+ H.c., (Dl) 

r,j ,s, s' 

where the phases 0j are the ones defined in Eq. © and 
we introduced a generic 2 x 2 unitary tunneling operator 
of the form O ss i = e laa . Here a is a generic position- 
independent vector and a is the vector of the Pauli matri¬ 
ces. In this case the canonical time-reversal transforma¬ 
tion must be defined as T = cr y K., where K, is the complex 
conjugation. It is easy to see that TOT ^ = O , therefore, 
once again, we obtain TH 2 ({9})T^ = H({—0}). Analo¬ 
gously to the spinless case, Eq. © holds. Hence, also for 
H 2l the physical time-reversal symmetry reads T = U y T, 
but, differently from the spinless case, we have T 2 = —1. 

Concerning the inversion symmetry, beside mapping 
r into —r, the parity must flip aa into —aa. There¬ 
fore P = V = i(3aV where (3 is a unit vector orthogo¬ 
nal to a and V describes the coordinate transformation 
© corresponding to V in the spinless case. In this way 
VOV t = O t such that VH 2 V t = H 2 . Also in this case 
[T,P} = 0. 

The Hamiltonian ( |D1| ) describes a Weyl semimetal 
with eight Weyl points in the Brillouin zone. With re¬ 
spect to the spinless case, the Weyl points are displaced in 
momentum space by the two vectors ± (|a|, |a|, |a|), such 

that k 0 = (±7r/2 + s|a|, ± 7 r /2 + s|a|, 7 r/ 2 -(-s|a|), 

where s = ±1 describes the two eigenstates of aa. For 





















|a| 7 ^ 0 , 7 r /2 this implies that the Weyl points do not 
overlap in momentum space and Eq. (D1) defines a spin 
1/2 PT-invariant Weyl semimetal with T 2 = — 1. 


The transformation of the Weyl points under PT can 
be easily obtained by generalizing the spinless case and 
considering that P maps s into — s. 


* for correspondence: michele.burrello@mpq.mpg.de 
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